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1 Introduction 

The study of the behavior of sohitions of ODEs often benefits from deciding on a convenient choice 
of coordinates. This choice of coordinates may be used to "simplify" the functional expressions that 
appear in the vector field in order that the essential features of the flow of the ODE near a critical 
point become more evident. In the case of the analysis of an ordinary differential equation in the 
neighborhood of an equilibrium point, this naturally leads to the consideration of the possibility 
to remove the maximum number of terms in the Taylor expansion of the vector field up to a 
given order. This idea was introduced by H. Poincare in [25j and the "simplified" system is called 
normal form. There have been several applications of the method of normal forms particularly in 
the context of bifurcation theory where one combines between the method of normal forms and 
the center manifold theorem in order to classify bifurcations [9]. This approach was extended to 
control systems in continuous-time by Kang and Krener {\17\, see also [18j for a survey) and Tall 
and Respondek ([23], see [23] for a survey), and by Barbot et al. [2] and Hamzi et al. in discrete- 
time [15^116]. The center manifold theorem was extended to control systems by Hamzi et al. |10|,lllj 
and combined with the normal forms approach to analyze and stabilize systems with bifurcations 
in continuous and discrete-time IHl (Hj . 

On another side, even though in many textbook treatments (see eg [9]) the emphasis is on the 
reduction of the number of monomials in the Taylor expansion, one of the main reasons for the 
success of normal forms lies in the fact that it allows to analyze a dynamical system based on 
a simpler form and a simpler form doesn't necessarily mean to remove the maximum number of 
terms in the Taylor series expansion. This observation, led to introduce the so-called "inner-product 
normal forms" in [HJ [T9| [7]. They are based on properly choosing an inner product that allows to 
simplify the computations. This inner-product will characterize the space overwhich one performs 
the Taylor series expansion. The elements in this space are the ones that characterize the normal 
form. Our goal in this paper is to generalize such an approach to control systems. 

In section §2, we review some results about normal forms. In section §3, we develop a new 
method for deriving normal forms for control systems. 



'Parts of this work were done while at Department of Mathematics, Duke University, Durham, NC 27708, USA. 
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2 Normal forms near equilibria of ODEs 



In this section we briefly review some results on normal forms near equilibria of nonlinear ODEs. 
Consider the nonlinear ODE in IR" 

x = Ax + fix), (2.1) 

with / G C"'"'"^(IR'^;IR'^), /(O) = and A = §|U=o is in real or complex Jordan form. Without loss 
of generality the latter condition can be met by application of a linear coordinate transformation. 
The goal is to find a change of coordinates 

x = m, (2.2) 

with G C"'(IR"; IR"") in a neighborhood of the origin, such that the Taylor expansion of (12. ip is 
simple, making essential features of the flow of (j2.ip near the equilibrium x = more evident. The 
desired simplification of (12. Ih will be obtained, up to terms of a specified order, by constructing a 
near identity coordinate transformation from a sequence of compositions of coordinate transforma- 
tions of the form (j2.2p with 

ay) = exp(eW)(y) = y + ^^(y) + 0{\y\'^'), (2.3) 

where y G IR" is close to zero, ^['^J e {k >2), the vector space of homogeneous polynomials 
of degree A; in n variables with values in IR", and exp{^^^^) denotes the time-one flow of the ODE 
y = ^^'^'(y)- We consider a formal power series expansion of / in (12. Ih and write 

/(x) = /Pl(x) + /[3l(x) + ..., (2.4) 

with /[''I G H^. From ([231) we obtain 

r\y) = y-^^'Ky) + Oi\y\'"'). (2.5) 

Substituting ([22]), ([23]) and ([23]) in (HT]), we get 

y = Ay + --- + /['^-il(y) + /W(y) - (LAeW)(y) + 0(|y|'^+^), (2.6) 
with the Lie derivative defined on vector fields / as 

(LAf){y):=^Ay-Afiy). (2.7) 

In the present context La is also known as the homological operator. 

The Lie derivative leaves H!^ invariant, : HI^ H^. We denote its range in as TZ^ and 
let denote a complement of TZ^ in 

H^ = n''®C\ k>2. (2.8) 

We define a normal form of / of order r as a Taylor expansion of the vector field with linear part 
and terms f^^^ G for 2 < k < r. We may associate the choice of complement to an inner 
product on H!^, for which it is the orthogonal complement of IZ'' in H!^, i.e. C'' := (TZ^)'^. 
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A convenient choice of inner product was introduced by Belitskii [3] , Meyer [19] and Elphick et 
al. , enabling the characterization of expression of as the kernel of the Lie derivative of A* (the 
adjoint of linear part A of the vector field at the equilibrium). Denoting monomials in shorthand 
notation as ^" with := ill ■ ■ ■ in^. , we define an inner product on polynomials 

p{^) = '^Pex^i q{x) = y^qmx"^, as {p,q) = '^m\pmqm- (2.9) 

Cm m 

For vector polynomials we define the corresponding inner product as the sum of the inner products 
between the polynomials of corresponding vector components. The inner product ()2.9p with T G 
gl{n,Wi) and T* denoting its adjoint (with respect to the standard inner product on IR") satisfies 

mm 

{poT,q) = {p,qoT*). (2.10) 

Accordingly, one obtains that the adjoint of La on HI^ with the above defined inner product satisfies 
the following relation [3l [7] 

{LaT = La*. (2.11) 

By application of the Fredholm alternative, it follows that iJZ^)-^ = ]iei{L*^\uk). In combination 
with (12. lip , this leads us to 

C'^ = ker(L^*|^.), 

as a result of which nonlinear elements of the normal form g satisfy the linear PDE 

LA^g = 0. (2.12) 

This PDE can be solved explicitly using the method of characteristics (for more details on this 
method, see for example [6]). 

We recall that since La* is a Lie derivative, it follows that the nonlinear elements of the normal 
form commute with the group 

G = {exp(yl*t) I t € M}. (2.13) 

We finally note that 

ker(L^.) = ker(LX.) n ker(L^. ), (2.14) 

where A = Aa + An is the Jordan-Chevalley decomposition of A in its (mutually commuting) semi- 
simple and nilpotent parts. As A* commutes with ^4^ but not with An (if nonzero), if A is not 
semi-simple, only a subgroup of G (as defined above) is a symmetry group of the normal form. In 
general, with the above choices made, the normal form is equivariant with respect to the group 



Gs = {exp(A*t) I t E IR}. (2.15) 
The appearance of this symmetry group is an important feature. 

3 Normal Forms of Nonlinear Control Systems 

The object of this section is to extend the normal form theory set out above to nonlinear control 
systems. We consider the nonlinear control system 

X = fix), (3.16) 
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with X = {x, u)'^ G IR" X ]R"^ and u € IR'" representing the control. In lowest linear order Taylor 
expansion, the control system takes the form 

X = Ax + 0(|xp), 

with A:={A 5 ) and A := ^U=o, B := ^k=o • 

We consider the effect of coordinate transformations of the form 

x = p{y) = exp(pW)(y) = y + + 0{\y\''+'), (3.17) 

where y = {y, v)'^ and € n with 

^n,m ■~ {(.Px^Pu^) I Px^ S -^n' pL ' ^ B^n+m,in}' (3.18) 

with H^j^^ ^ denoting the vector space of homogeneous polynomials of degree k from IR""''™ to 

IR'". The skew product form of p^^^ G S^,mi P^^\v) = ijP^x\y),P^u\y))'^ , guarantees that the control 
system is transformed to another control system of the same type. If would depend on u then 
the coordinate transformation would introduce a relationship involving n. We obtain 

y = Ay+... + f^^-^\y) + /W(y) - {C^p^^^){y) + 0(|y|'^+i), (3.19) 

where the homological operator : „ — >■ -ff^+m „ has the form 

{C-AP^^^m = Dp^^\y)Ay - ^pW(y) = (LApL")(y) + Dp^x\y)Bn - Bp^^\y). (3.20) 

We recognize in this expression the Lie derivative La, that is equal to Ca in case i? = 0. Indeed, 
f{y,0) (the part of / that does not depend on u) can be put into a G^-equivariant normal form, 
using coordinate transformations of the form p{y) = {ex.p{px){y),u) only. 

We now proceed to characterize a normal form by a (choice of) complement of the range of 
Ca- In order to do so in analogy to the theory developed for ODEs, we temporarily take the 
viewpoint as if the coordinate transformation would be for the ODE (x,n) = {f{x,u),h{x,u)), for 
some h : ST' with L>^/i(0,0) = and L'„/i(0,0) = 0. The homological operator for the 

latter ODE, with coordinate transformations of the form p.l7p takes precisely the form of the Lie 
derivative Laq, with 

A f ^ B \ 
- V j ' 

so that Laq = {Ca,0) and Ca •= '^^Aq, with vr : denoting the canonical projection 

7r(x, u) := X. 

We may thus choose the complement of the range of Ca as the projection under vr of the 
orthogonal complement C'^ to the range of Laq taken with respect to the inner product (j2.9p with 

& := {q G I {q, La,p) = 0, Vp G S^^}. (3.21) 

By the Predholm alternative we have 

{q,LA,p) = {LA'q,p), (3.22) 
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so that this complement takes the form 

:= {q G I ^ «^)^}. (3.23) 

By the definition of the inner product (12. 9p , 

(^n, = e I g(x,o) = o}, 

i.e. the subset of vector polynomials in -ff^+n which each constituting monomial contains 
a factor lij, i = 1, . . . ,m. The complement to the range of C_a characterising the corresponding 
normal form is vrC'^. By writing out the relevant operators, the following result follows immediately. 

Theorem 3.1 (Control normal form) Consider a finite order in Taylor expansion of the vector 
field defining the control system /i3.16\) . 

N 

/(x) =Ai + Y. /"^(^) + 0(1^1'+'), wzth /W G 

k=2 

By a choice of coordinates, the nonlinear parts /t'^l can be made to satisfy 

^^•/t"(^,0) = (3.24) 

where 

£^./W(x) := Z),/W(i)-4*x - A*f^'^\x), (3.25) 

and X = (x, u). 

Remark. We note that by restricting first to coordinate transformations that do not involve 
we can achieve G^-equivariance of the control system to any desired order. Then we can refine the 
normalization further using G^-equivariant coordinate transformations that preserve this equivari- 
ance. <i 



4 Illustrations 

4.1 Linearly Controllable Case 

To illustrate this method, consider the nonlinear control system S in (|3.16p with one input, i.e. 
m = 1, and assume that its linearization is controllable. From linear control theory we know that 
there exists a linear change of coordinates and feedback that allows to transform the linear part in 
the Brunovsky form, i.e. 



A 



/ 1 

1 



V 



\ 


1 

/ 






V 1 / 



(4.26) 
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In this case, the PDE (j3.25p becomes 

dpi 



__ dpi_ dpi 

0X2 OXn OU 

dp2 , , 5p2 , dp2 

Xl^ \ h Xn-l^ h Xn-^ Pi 

0X2 OXn dU 



dpn , , dpn , dpn 

Xit; \ hXn— Pn~l 

0X2 OXn OU 



(4.27) 



that we'll solve using the method of characteristics. 

Theorem 4.1 Consider the nonlinear control system S given by i3.16\) . There exist a change of 
coordinates and feedback {3.11^ such that S writes as 

xi = X2 + ^i{ii, - ■ ■ Jr+l] 



X2 



X3 + '^2{h,- 



,ir+l) + I Plix,u) — , 
J Xi 



(4.28) 



U + ^n{h,--- Jr+l)+ J Pn 



. dX2 

-i{x,u) , 

Xl 



with ^i{ii 



,£n) are functions satisfying 



for p = 0, ■ ■ ■ ,n — i, 



xi=0 



(4.29) 
(4.30) 



n—p 



X 1 ^— ^ 

and li{x) = Xl, hix) = - X1X3, • • •, li{x) = -xj + 2_^{-^)^Xi-kXi+k for i = 2, ■ ■ ■ ,r + 1, and 



k=l 



r is such that r = n/2 if n is even, and r = {n — 1) /2 if n is odd (here, xq = and x. 



n+l 



U) 



Proof 

In the n— dimensional space of the variables xi, X2, - ■ - jXn we determine the curves Xj = Xj(s) 
in terms of a parameter s by means of the system of ordinary differential equations that represent 
the characteristic curves 

dxi 
ds 

dX2 

ds 





Xl 



(4.31) 



dx. 
~d, 



is 
du 

ds 



Xn~l 
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Along the characteristic curves and using the chain rule, the systems of PDEs (j3.25p writes as 

dpi dxi dpi ^ dx2 dpi ^ ^ dxn dpi ^ du dpi ^ 

ds ds dxi ds dx2 ds dxn ds du 

dp2 dxi dpi dx2 dpi dxn dpi du dpi 

ds ds dxi ds dx2 ds dxn ds du (4.32) 

dpn dxi dpi dx2 dpi dxn dpi du dpi 

ds ds dxi ds dx2 ds dxn ds du ^" ^ 

Hence, along the characteristic curves defined by (j4.3ip . the systems of PDEs (|3.25p transforms 
into a set of ODEs 



dp 



d% 



1 

'P2 







ds (4.33) 



dpn 

ds 

This system of ODEs can be solved explicitly 



Pn~l 



Piis) = Cl 

P2{s) = + y Pi{s)ds 



Pn{s) = Cn + j Pn-l{s)ds 

The "constants of integration", Cj, are the constants along the characteristic curves which are 
the trivial first integrals of the system (j4.3ip . One can check that they are given by ^i(x) = xi, 

x^ 1 
^2(3;) = — — 2:1X3, • • •, ti{x) = -xf + 2~]{~^)^ ^i-kXi+k for i = 2, ■ ■ ■ , r + 1, and r is such that 

k=l 

r = n/2 if n is even, and r = (n — l)/2 if n is odd (for notation convenience, xq = and Xn+i = u). 
From ()4.3ip , we hav^ll ds = ^ , and the solution of ()4.33p is given by 

Pi{x,u) = ^>i(^i, • • • ,4+1) 

P2ix,u) = ^2{il,---,ir+l) + Jpi{x,u)^ 

(4.34) 

Pn{x,u) = <^n{h,---Jr+l) + fPn-l{x,u)^ 

where $i(^i, • • • ,in), ^ = 1; ■ ' ' are functions of the variables ii, - ■ ■ ,in and are thus constants 
along the characteristic curves define in (|4.3ip . Since ^i{£i, - ■ ■ ,£n) and qi{x,u) satisfy the condi- 
tions 



= for p = 0, ■ ■ ■ ,n 

xi=0 



^We can also use ds = — and in this case the normal form will be parametrized by Xi+i. We can also 
parameterize each component with a different parameterization. 
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thus pi, ■ ■ ■ ,Pn-i are divisible by xi. Hence pi{x,u), ■ ■ in (j4.34p are polynomials. 



4.1.1 Example 

Consider a two dimensional system with controllable linearization. In this case, the linear part of 
(|3.16p writes as 

xi = X2 (4.35) 
X2 = u (4.36) 



and the PDE ([3:25]) writes as 



dpi dpi 

Xl- hX2^ = 

0x2 ou 
dp2 , dp2 

0x2 ou 



Hence, we get 



dpi 
ds 

dp2 
ds 



with 



(4.37) 



: 0, (4.38a) 
Pi = (4.38b) 



(4.39a) 

Xl (4.39b) 

X2 (4.39c) 

(4.39d) 



dxi 
ds 

dX2 

ds 
du 
ds 

dpi 
ds 

^-Pi = (4.39e) 
ds 

We thus deduce the following parametrization of the solution 

Xl = xifi (4.40a) 
X2 = xifiS + X2,o (4.40b) 
= + X2,0S + X3,o (4.40c) 

The first integrals are £i{x,u) = xi and i2ix,u) = 2xiu - x\. From (I4.39dl) - (l4.39e|) we deduce 
that 

Pi{x,u) = ^i{xi,2xiu-x\) (4.41) 
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P2{x,u) = j pi{t)dt + ^2{xi,2xiu - xl) (4.42) 

We can use either ()4.39bp or (|4.39c|) to express the normal form as a function of X2 or u. For 
example, using (j4.39bj) we deduce that dt = — -. Moreover, using ()4.29p . we obtain conditions on 
^i{h,(.2) and i = 1,2, 



At the quadratic level these conditions imply that 

Mh,i2)\x^=o = ^nxl + 0{x,uf (4.43a) 

^2{ii,i2)\x,=o = ki4 + 4>i2i2x,u-xl) + 0{x,uf (4.43b) 

Hence 

Pi{x,u) = (t)iix\ + 0{x,uf (4.44a) 

P2{x,u) = cl)iiXiX2 + (t)iixl + 4>i2{2xiu - x1) + 0{x,uf (4.44b) 

Hence the normal form has the form 

xi = X2 + (t)iix\ + 0{x,uf (4.45a) 
X2 = u + (j)iiXiX2 + (t>iix\ + (j)i2{'2.xiu - x\) + 0{x,uf (4.45b) 

4.2 Systems with Uncontrollable Linearization 

Now, consider the nonlinear control system S in (j3.16p with one input, i.e. m = 1, and assume 
that the system has r uncontrollable modes. From linear control theory we know that there exists 
a linear change of coordinates and feedback that allows to write the linear part as 

z = Aiz + 0{z,x,uf, (4.46) 
X = A2X + B2U + 0{z,x,uf (4.47) 

where z € W^, x G IR("-^)^i, Ai G IR'^^^ and (^2,^2) S IrC'^-O ><(«-'■) x IR("-^)xi are in the 
Brunovsky form. 

In this case, ^=(^^ ^ ),a; = {z,x,u)'^ in the PDE ()3.25p . Let's note that when r = 



^ A2 B2 ^ 

we recover the case in the preceding section and we can find a general explicit solution. However, 
when r 7^ a general solution is not as easily found and depends on Aq. We'll illustrate the method 
through an example. 
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4.2.1 Example 



Consider the system whose hnear part writes as 



z = 0(2:, 
xi = a;2 + 0(2:, x, n)^ 



(4.48) 



This system has uncontrohable hnearization and the uncontrollable dynamics corresponds to the 
z— dynamics. 

The elements of the normal form satisfy the PDE 



^^2^— 

0x2 ou 

dp2 , dp2 

^7, ^ ^2^— 

0x2 ou 

dp3 dp3 



dx-^ 



+ X2 



du 



The equation of the characteristics is 



dpi 

ds 
dp2_ 

'P3 

ds 



-P2 




P2 









(4.49) 



(4.50) 



The characteristic equations are 2; = ci, xi = C2, X2 = cis + C3, u 
parametrize by X2 or u by writing ds = or ds = 
The solution of the system of PDEs (j4.49p is 



+ C3S + C4. We can either 



Pi 

P2 
P3 

The normal form is thus given by 



'^o{z,xi,xl 



2zu) 



^'i(z, xi,X2 — 2zu) 
■^2{z,xi,xl-2zu) + fp2{z,x,u)^ 



(4.51) 



z 
±1 

X2 



2zu) 



'^o{z,xi,xl 
X2 + ^'i(2;,xi,X2 — 2zu) 
u + ^'2(2, xi,xl- 2zu) + / P2{z, X, u)^ 



(4.52) 



5 Concluding remark and future extensions: 

Given the preceding, one could think about hyper normal forms where instead of normalizing with 
respect to the linear term, one normalizes the quadratic term with respect to the linear term, then 
normalize the cubic term with respect to the sum of the linear and quadratic terms, and so forth. 
This direction has been fruitful for systems without control \2\\ [22] and its extension to the control 
case is the object of future research. Several other extensions are possible for this work. One could 
think about characterizing completely the normal form in the case of systems with uncontrollable 
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linearization, developing the Hamiltonian case, and computing the coefficients in the normal form 
directly from the original system. 

Acknowledgement: The first author is thankful to Prof. Murdock for useful comments and 
for pointing out to references [H [22]. BH also thanks the European Union for financial support 
received through an International Incoming Marie Curie Fellowship. 
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